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Abstract 


The  paper  contains  an  analytical  solution  to  the  title  problem  by  means  of  the  semi¬ 
membrane  theory  of  cylindrical  shells.  Results  include  the  crack  opening  displacements  as  a 
function  of  circumferential  angle,  the  fracture  parameter  J,  and  the  effects  of  the  crack  on  the 
compliance  of  the  pipe. 


INTRODUCTION 

Ideally  the  failure  load  of  a  cracked  pipe  could  be  calculated  from  theory  and  the  results  of  a 
small  number  of  tests  on  simple  specimens.  Unfortunately  this  is  not  the  case  nor  is  it  likely  to  be 
any  time  in  the  near  future.  A  good  idea  of  the  "state  of  the  art"  can  be  gained  from  reading  recent 
reports  on  research  projects  conducted  by  General  Electric  [1],  Battelle  [2],  and  the  Proceedings  of 
an  A.S.M.E.  Symposium  [3].  The  fundamental  theory  of  nonlinear  fracture  mechanics  is  still 
being  developed.  How  it  applies  to  real  structures,  pipes  in  particular,  is  not  a  routine  matter. 
Engineering  procedures  based  on  available  knowledge  plus  a  considerable  amount  of  expensive 
testing  must  be  used  to  get  information  for  maintenance  or  design  purposes. 

Frequently  the  practice  is  to  apply  experimentally  determined  factors  to  results  for  cracks  in 
flat  plates  to  correct  for  curvature  effects  or  to  make  adjustments  to  results  from  elastic  solutions 
for  cracks  in  cylindrical  shells  to  account  for  plasticity.  Such  procedures  are  not  unreasonable  but 
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leave  unanswered  such  questions  as:  can  tests  on  small  diameter  pipes  be  trusted  to  predict  what 
w  ill  happen  to  large  diameter  pipes?  It  seems  that  more  effective  use  could  be  made  of  cylindrical 
shell  theory  without  resort  to  finite  element  methods.  The  present  paper  is  directed  toward  that  end 
in  the  case  of  circumferential  through-cracks. 

The  analysis  to  be  presented  has  some  limitations  of  its  own,  however.  It  applies  when  the 
pipe  material  is  such  that  the  Dugdale  model  for  plasticity  effects  is  appropriate.  It  applies  only 
when  the  crack  is  sufficiently  long  for  the  semi-membrane  theory  of  cylindrical  shells  to  be  valid. 
For  short  cracks  with  a  Dugdale  zone  shallow  shell  theory  should  be  used,  and  that,  strangely 
enough,  is  a  more  difficult  problem.  It  could  be  solved  but  no  one  has  yet  done  it.  Most  interest 
seems  to  be  in  cracks  which  extend  at  least  an  eighth  of  the  circumference  and  that  is  long  enough 
for  the  present  analysis  to  be  valid.  For  some  further  discussion  see  [4].  The  analysis  here  is 
limited  to  the  bending  case  because  that  seems  to  be  most  important  technically.  Tension,  torsion, 
or  combined  loading  cases  can  be  done  by  the  same  methods. 

Results  of  the  sort  obtained  here  could  be  obtained  by  finite  element  methods  and  without 
need  for  the  Dugdale  model  but,  of  course,  at  the  expense  of  much  more  computing.  The  principal 
advantage  of  the  present  approach  is  that  the  results  come  out  in  closed  form  in  terms  of  elementary 
functions.  Some  numerical  root-finding  is  required,  but  programming  to  do  parameter  studies 
would  be  trivial.  The  analysis  includes  linear  elastic  fracture  and  plastic  collapse  as  limiting  cases. 
A  theory  of  tearing  based  on  the  present  analysis  and  an  additional  physical  hypothesis  has  been 
constructed  and  will  be  the  subject  of  a  subsequent  paper. 

FUNDAMENTAL  EQUATIONS 

Dimensionless  displacements,  stress  functions,  and  stress  measures  are  related  to 
corresponding  dimensional  quantities  (with  an  overbar)  by  the  relations 
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(u,  v,  w)  = 


oFR 


(u,  v,  e_2w) 


(Xz>Xfl’  V)=  cFR2h£2(Rx  ,  Rxe »  e~2V) 


(Nz>Ne,Nz9)=oFh(Nz,Ne,N2e)  (1) 

(Mz,  Mg.  Mze )  =  (5pRhe2(  Mz,  Mg,  Mzg ) 

M=n  R2haFCT 

Here  aF  is  the  flow  stress  of  the  material,  h  is  the  pipe  thickness,  R  is  the  pipe  radius,  and  E  is 
Young's  modulus,  M  is  the  applied  bending  moment  on  the  pipe,  a  is  a  dimensionless  load 

parameter  and  e  is  a  small  parameter  given  by 

e2=  (h/R)[12(l— v2)]"1/2  (2) 


The  analysis  in  the  present  paper  is  similar  to  that  in  [4]  where  it  was  found  that  the  most 
significant  results  are  obtained  by  assuming  a  semi-membrane  state  of  stress  in  the  cylinder  and 
neglecting  the  small  contribution  of  edge  effects.  The  solution  thus  obtained  is  correct  except  in 
the  immediate  vicinity  of  singularities  such  as  the  crack  tips.  Nevertheless  quantities  such  as  stress 
intensity  factors,  overall  compliance,  and  so  forth,  important  to  fracture  mechanics,  turn  out  to  be 
correct.  A  semi-membrane  state  is  derived  from  two  complex- valued  functions  O  and  <p  each  of 
which  satisfy  the  characteristic  equation 

(0+<X>)"-iE-20”  =0  (3) 

Here  primes  and  dots  refer  to  differentiation  with  respect  to  z  or  0  respectively  where  zR  is  the 
axial  distance  measured  from  the  cracked  section  and  0  is  the  circumferential  angle  measured 

from  the  middle  of  the  crack.  The  functions  O  and  9  are  related  by 

<J>"  =  e2cp  (4) 

Dimensionless  displacements  and  stress  functions  are  given  in  terms  of  C>  by  the  formulas 

e2u  =  O'  £2XZ  = 

e2v=-0  £2X0=i&  (5) 


w  =  -<J>+  i9 


y  =  i0>+9 
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Dimensionless  membrane  and  bending  stress  measures  are  given  in  terms  of  9  by 


Nz  =9 

Mz  =  -i(9"  +  v9) 

Nq  =  <P" 

M0  =  -i  (9+  v9") 

(6) 
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In  these  formulas  the  physical  variables,  displacements,  etc.,  are  to  be  understood  as  the  real  part 
of  the  expressions  on  the  right-hand-side  of  the  equation. 

ANALYSIS 

Because  of  symmetry  boundary  conditions  can  be  applied  on  the  section  z  =  0  .  In 
physical  terms  the  crack  edges  are  free  of  stress.  A  Dugdale  zone  is  assumed  to  exist  beyond  the 
crack  tips  for  an  as  yet  undetermined  distance.  Within  the  Dugdale  zone  the  axial  membrane  stress 
equals  the  flow  stress.  The  axial  bending  stress  is  ignored  on  the  basis  of  the  assumption  that  the 
axial  stress  equals  the  flow  stress  throughout  the  thickness  of  the  pipe.  For  a  sufficiently  large 
applied  bending  moment  compressive  stresses  on  the  side  of  the  pipe  far  from  the  crack  could 
exceed  the  flow  stress  and  so  a  compressive  Dugdale  zone  (of  unknown  extent)  is  admitted  to 
prevent  it.  All  remaining  boundary  conditions  on  z  =  0  follow  directly  from  the  symmetry  of  the 
problem.  In  the  analysis  the  pipe  is  assumed  to  be  infinitely  long  and  subjected  to  a  bending 
moment  at  infinity  in  such  a  way  that  the  crack  is  symmetrically  located  on  the  tension  side. 
According  to  the  semi-membrane  theory  only  two  boundary  conditions  apply  on  any  segment  of 
the  section  z  =  0  .  In  other  words  two  of  the  possible  four  conditions  are  redundant.  Boundary 
conditions  on  stress  measures  are  more  conveniently  expressed  in  terms  of  conditions  on  stress 
functions  obtained  in  accordance  to  procedures  in  [5].  Because  of  symmetry  the  boundary 
conditions  need  only  be  stated  in  the  range  0  <  9  <  7t  .  In  terms  of  dimensionless  quantities  they 
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e2xe  =  0  (O<0<a) 

=  sin(0 -a)  -  (0-a)  (a<0<P)  (8) 

=  sin  0  —  (7T  -0)  +  -y  acos0-  Gj  sin  0  (y<  0  <7t) 

u  =  0  (p  <  0  <  Y)  (9) 

Gj  is  a  constant  of  integration.  For  positive  values  of  0  the  crack  extends  from  0  =  0  to  0  =  a  , 
the  tension  Dugdale  zone  lies  in  a  <  0  <  P  ,  the  compressive  Dugdale  zone  (if  any)  lies  in 
Y  <  0  <  tc  ,  and  similarly  for  negative  values  of  0  .  By  means  of  equation  (5)  the  boundary 
conditions  (on  z  =  0)  can  be  restated  in  terms  of  $ 


R  {i<J>c}  =  0 

an  0 

O'c  =  0 

P<0  <  Y 

R  (iOc}  =  0 

0^0  <a 

=  1-  ^ (0-a)2- cos (0-ot) 

a<0  <  P 

=-  (l+cos0)+  -i (re  — 0)2  +^-asin0  +  Gj  cos0  +  G2 

Y <0  -it 

G2  is  a  constant  of  integration.  The  subscript  c  has  been  introduced  to  indicate  that  these 

conditions  apply  to  the  complete  solution.  The  complete  solution  C>c  is  a  superposition  of  a 

function  O  which  can  be  assumed  to  vanish  at  infinity,  a  beam  solution  (which  accounts  for 

bending  by  moments  at  infinity)  given  by 

<pg  =  -i[e2z2cos0  +  i(2cos0-0sin0]o  (11) 

and  a  collection  of  four  elementary  solutions  given  by 

<i>E=  ia-ebz+ ic  cos0  +  edz  cos0  (12) 

which  in  physical  terms  mean  a  null  solution  (does  nothing),  a  real  axial  translation,  an  imaginary 
lateral  translation,  and  a  real  rotation.  The  constants  a  ,  b  ,  c  ,  and  d  are  real,  and  as  yet 
undetermined.  In  case  there  is  no  compressive  Dugdale  zone  then  y=n  and  the  third  of  equations 

(8)  becomes  inoperative.  The  boundary  conditions  for  the  reduced  problem  read 

fl{iO'}  =  0  all  0  (13) 
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o'  =  e(b-d  cos0) 

P  <  0  <  Y 

(14) 

R  { iO}  =  I(2cos0 -0  sin0)a+  a  +  c  cos0 

0  <  0  <  a 

=  i(2cos0  -0  sin0  )c+  a  +  c  cos0 

+  l-^(0-oc)2-cosf0-a) 

a  <0  <  p 

(15) 

=  4l2cos0-(0-Jt)sin0]a+  a+  ccos0 

-(1+COS0)  +  (Jt  -0)2+  G!  COS0+  G2 

y  <0<  it 

According  to  [6]  the  necessary  and  sufficient  conditions  for  O 

to  vanish  at  infinity 

are 

rn  ••  r*  •• 

Od0=O  .  O  cos 0  d0  =  0 

Jo  J  0 

(16) 

and  in  such  a  case  O'  is  accurately  approximated  by 

6'  =  -i3/2£(o+ior 

(17) 

Now  put 

O  =  i1/2H(0)  ,  on  z  =  0 

(18) 

then  by  (17) 

O'  =  e(H+^H)  ,  on  z  =  0 

(19) 

Equation  (13)  will  be  satisfied  if  H  is  real.  Equation  (14)  becomes  a  simple  differential  equation 

for  H  in  the  range  P  <  0  <  y  ,  while  (15)  gives  H  explicitly  over  the  rest  of  the  range  of  0  . 

Complete  results  for  H  are  as  follows 

H=  n/2  (a  +  ccos0)- V2  (cos0-^0sin0)a 

H  =  -V2  (a  +  c  cos0)-  V2(cos0  -  J  0  sin  0)a 

O<0  <a 

(20) 

-VI  (1-^  (0-a)2-cos (0-o)] 

a  <0<p 

(21) 

H=  2(b  +  dcos0)+ 4Pcos  x+ 4Q  sin  x 

P  <  0  <  y 

(22) 
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H  =  -V2[a+G2+(c+Gi  )cos0] 

-a/2[cos0+  ~  (7t-0)sin0lc 

+  V2[1+cos0-|(t:-0)2]  J<Q<n  (23) 

On  physical  grounds  the  displacements  and  stress  functions  must  be  continuous  at  a  ,  p  , 
and  y .  This  means  that  H,H,  and  H  must  be  continuous  at  p  and  y.  Continuity  at  a  is 
automatic.  The  extent  of  the  Dugdale  zones  is  determined  by  the  condition  that  xz  (°r  H)  be 
continuous  at  p  and  y .  Recall  that  in  plane  stress  applications  of  the  idea  an  analogous  condition 
must  be  imposed  to  avoid  infinite  stresses.  In  the  present  context  a  discontinuity  in  H  means 
(mathematically)  the  presence  of  a  concentrated  load  but  it  should  be  understood  that  semi¬ 
membrane  theory  is  an  asymptotic  approximation  in  which  a  large  stress  acting  over  a  short 
distance  gets  replaced  by  a  concentrated  load. 

To  simplify  the  discussion  a  bit  consider  the  case  in  which  there  is  no  compressive  Dugdale 

zone.  Then  y=  7t  and  Q  =  0  .  There  are  six  unknown  constants  namely,  a  ,  b  ,  c  ,  d  ,  P  ,  and  p. 

There  are  six  equations  namely,  the  four  continuity  conditions  on  H  and  its  first  three  derivatives 

and  the  two  conditions  (16)  which  can  be  written  in  the  equivalent  form 

rn  .. 

(H+H)d0  =  0  ,  Hcos0d0  =  O  (24) 

Jo  Jo 

As  it  turns  out  the  continuity  conditions  involve  a  ,  b  ,  c  ,  and  d  only  in  the  combinations 

A=a+V2b  ,  C=c  +  V2d  (25) 

As  an  algebraic  system  the  four  continuity  conditions  are  four  linear  equations  for  the  three 
constants  A  ,  C  ,  and  P .  The  condition  that  the  equations  be  compatible  is 

c  =  {[sin  P- sin  a  +  (p-a  )cos  p] sin  y  +  V2  ((5 -a  )sin  Pcos  y} 

+  ["^(P+  3  sin  P  cos  P)sin  y  +  V2  sin2  Pcos  yl  (26) 


where 
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y  =  (n  —  p)/V2 

This  is  a  transcendental  equation  for  the  determination  of  P  given  the  load  parameter  a  and  the 
crack  length  a  .  Only  A  ,  C  ,  and  P  are  required  to  determine  the  crack  opening  displacements 
uc(O,0) ,  and  (according  to  the  Dugdale  model)  the  fracture  parameter  J  =  2opuc(O,0) .  However, 
b  and  d  are  needed  to  calculate  the  extension  and  rotation  of  the  pipe  (over  and  above  beam 
theory)  at  sections  remote  from  the  crack.  These  deformations  relate  to  the  increased  flexibility  of 
the  pipe  due  to  the  presence  of  the  crack.  The  constants  in  question  are  determined  by  use  of  (24). 
Results  for  A  ,  C  ,  P  ,  Q  ,  and  0  are  listed  here  for  the  case  y  <  K  .  If  no  compressive 


Dugdale  zone  is  present  replace  y  by  n  . 

Q  =  0  sin  y  -  Ot  -  y)  (27) 

Psin  x  =  Qcosx- asin  P  +  p  -  a  (28) 

x  =  (y_p)/V2  (29) 

A  =  acosP+  1  (P-a)2-  V2(Pcos  x+  Qsin  x)  (30) 

CsinP=-^(7sin  P+Pcosp)a+p-a  +  sin(p-a)  (31) 

N,  =  [sin  p -sin  a  +  (P-a) cos p]  sin  x  -  V2  [jt  - y-  (p-a)cos  x]sin  p  (32) 

D*  =  ^P+3sinpcosp)sin  x-V2(sin  y-sin  Pcos  x)sin  p  (33) 

N2  =  [sin  y-  (tc  -  y)cos  y]  sin  x  +  V2[(7t  -  y)cos  x  -  (p-a)]  sin  y  (34) 

D2  =  ^2  (sin  ycosx-sin  P)sin y- ^[3  sin  ycos  y-  (n  -  y)]  sin  x  (35) 

0=  Nj/Di  =  N2/D(2  (36) 


Incase  y  =k  the  formula  o  =  Ni/Dj  reduces  to  (26).  The  value  of  0  for  which  the  axial 
membrane  stress  at  9  =  7t  first  reaches  the  flow  stress  is  given  by 

0=  (V2  sin  y -P+ a  )/(V2  sin  y  -  sin  p)  (37) 

where  y  =  (it  -  P)/v2  . 


I 


At  the  plastic  collapse  load 


Y  =  (3  =  (a  +  7t)/2 


(38) 


and  the  collapse  value  of  the  load  parameter  is  given  by 

a,  =  4  (cos  ^  - 1  sin  a)/7t  (39) 

L  2  2 

For  y  <  7t  the  expressions  in  (36)  must  be  used  simultaneously  to  determine  p  and  y 
given  a  and  a  by  some  numerical  root-finding  process.  Small  values  of  n  -  y  are  especially 
troublesome  because  the  second  expression  for  o  in  (36)  becomes  indeterminate  at  y  =  Jt .  A 

polynomial  approximation  for  (sin  z)/z  given  in  [7]  may  be  found  helpful. 

Results  for  a  ,  b  ,  c  ,  d  ,  and  Gj  follow  from 

n  a  =  (rt  -  p)A+  [sin  P~sin  y-(7t  -  y)cos  y]c+  -g[(p-a)3-  (k  -  y)3+  3(rc  -  y)4] 

+  [  2  sin  x+  V2  (jt  -  y)]P+  2(l-cosx)Q  (40) 

b  =  V2(A-a)/2  (41) 

G!Cosy=  1  +cosy-CcosY-^2P-[2cosy  +  -|(7t-y)siny]a  (42) 

7tc  =  -(P~Y  +  sin  p  cos  P -sin  ycos  y)C  -  (7t  -  y-  sin  ycosy^j 

-  2a[p  +  ^sinp  cosp  +  ^Pcos  2p  +K-  y-^  sinycosy+|(7t  -  y)cos  2y] 

-  2[  sin  P  -  sin  a  -  ^  ( p  -  a  )cos  a  -  ^  cos  p  sin  ( P  -  a ) 

+  sin y- 1  (7t  -  y)+^  sinycos y]+  2V2P[2  sin  y  -  2  sin  P  cos  x-  V2  cos  p  sin  x] 

-  2^2 Qf  V2  cos  y+  2  sin  p  sin  x-V2  cosP  cos  x]  (43) 

d  =  V2  (C  -  c)/2  (44) 

A  general  expression  for  the  complete  axial  displacement  (dimensionless)  on  the  cracked 
section  is  given  by 

euc(0,9)  =  H  +  ^H -b  +  dcos9  (45) 

Displacements  on  the  crack  are  obtained  from  (45)  by  using  expressions  for  H  valid  in  the  range 
0  <  9  <  a  .  The  result  is 
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euc(0,9)  =  ^  V2  [(3cos0  -^0sin0)a-A  +  Ccos0]  (46) 

The  dimensionless  crack  opening  displacement  8  =  2uc(0,  a)  is  given  by 

eS  =  V2[(3  cos  a  -  a  sin  a)a-  A  +  C  cosa]  (47) 

According  to  the  Dugdale  model  the  fracture  parameter  J  equals  the  flow  stress  times  the 

crack  opening  displacement  In  dimensional  form  the  result  is 

J  =  (^R8/E  (48) 

Values  of  the  field  variables  (stress  measures,  displacements,  etc.)  could  be  calculated  for 
values  of  z  other  than  z  =  0  by  solving  the  characteristic  equation  (3)  for  O  with  boundary 
conditions  according  to  (18)  and  (20)  to  (23),  possibly  by  means  of  Fourier  series.  At  present  it 
seems  there  would  not  be  much  interest  in  such  calculations.  Disturbances  in  the  stress  distribution 
due  to  the  crack  decay  exponentially  as  z  increases,  otherwise  the  pipe  behaves  as  a  beam.  How 
rapidly  these  effects  decay  depends  mostly  on  the  value  of  £  .  The  smaller  e  is  the  slower  the 
decay.  For  e  =  .1  the  effects  have  mostly  died  out  three  or  four  diameters  away  from  the  cracked 
section.  The  crack  does,  however,  cause  a  rotation  and  an  axial  displacement  which  persist. 

Dimensionless  axial  displacements  approach 

£uc  =  eoz  cos  0  —  b  +  d  cos  0  as  z  — >  oo  (49) 

The  term  with  z  corresponds  to  the  pipe  bending  as  a  beam  fixed  at  z  =  0  .  The  additional 
physical  rotation  and  axial  displacement  of  cross  j?ctions  are  given  by 


£2  =  aFd/E£ 

(50) 

U  =  -  Op  Rb/E£ 

(51) 

Results  from  the  present  analysis  should  reduce  to  those  from  linear  elastic  fracture 
mechanics  in  the  limit  Op  -»  .  Since  J  must  remain  finite  (48)  shows  that  8  ->  0  in  the  same 

limit,  as  expected.  The  extent  of  the  Dugdale  zone  P-a  must  also  become  small.  Somewhat 
lengthy  calculations  show  that  for  small  P-a 


£8  =  V2(P  — a  )2 


(52) 


o  =  ({3-a)/|isin  a 


where 


1+V2_  a  +  a  cot 2  a  -  cot  a 
^  cot  [(ti  -  a)/V2]  +  V2  cot  a 


(54) 


and 


^£)  =  -  e_1  V2  (p-a ) 

owe=a 


(z  =  0) 


(55) 


These  equations  imply 


eS  =  V2  o2(i2  sin2  a 


(56) 


It  follows  then  from  (48)  that 


J  =  ^r~  Cp02|l2sin2a 


(57) 


but  from  ( 1 )  repeated  here 


M  =  7tR2haFc 


(58) 


it  follows  that 


J  =■ 


V2 


7t  2EeR3h2 


M2q2  sin2  a 


(59) 


This  agrees  with  the  result  in  [8]  if  account  is  taken  of  the  relation 

I  =  2h  RJ 

between  the  I  in  [8]  and  J.  Equations  (48),  (52),  and  (55)  lead  to 

,=i'r2EeR(ilL 


(60) 


(61) 


which  with  J  =  K2/E  yields 


W  -  7-l'4cl'2nl'2p(  1  due 

K“2  e  R  E(r3TV 


0=ot 


(62) 


This  last  equation  shows  that  according  to  semi-membrane  theory  in  the  elastic  case  the  opening 
angle  at  the  crack  tip  is  a  measure  of  the  stress  intensity  factor.  Why  this  should  be  so  must  remain 
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somewhat  obscure  for  the  present  because  an  analytical  solution  to  the  full  shell  equations  valid  in 
the  neighborhood  of  the  crack  tip  is  not  available.  The  angle  here  must  be  the  limit  of  some 
quantity  obtained  from  such  a  solution. 

CONCLUDING  REMARKS 

Rather  complete  results  have  been  obtained  for  the  problem  of  through  circumferential 
cracks  in  pipes  with  use  of  the  Dugdale  model  and  according  to  thin  shell  theory.  One  would 
expect  the  theory  to  be  more  accurate  for  low  strain  hardening  materials  and  for  thin  pipes.  The 
thinness  requirement  is  not  too  severe  as  far  as  shell  theory  is  concerned.  Shell  theory  is 
reasonably  good  for  values  of  R/h  as  small  as  ten  and  not  grossly  in  error  for  R/h  as  small  as 
five  or  six.  See  the  paper  by  LeFort  et  al.  in  [3].  Errors  involved  in  the  use  of  the  Dugdale  model 
as  applied  in  the  present  paper  are  more  difficult  to  assess.  The  shape  of  the  plastic  zone  near  the 
tip  of  a  circumferential  crack  is  probably  quite  different  from  what  it  is  in  a  flat  plate.  Stresses  tend 
to  decay  more  slowly  in  the  axial  direction  and  equalize  more  rapidly  in  the  circumferential 
direction  than  they  do  in  a  plate.  One  might  be  tempted  to  argue  either  way  on  the  matter  but 
perhaps  it  would  be  better  to  await  the  results  of  more  refined  analyses  and/or  the  results  of 
experiment  A  sort  of  adaptation  of  the  line-spring  model  might  allow  the  theory  to  deal  with  strain 
hardening  in  a  more  realistic  fashion  but  there  seems  not  to  be  much  point  in  pushing  the  model  too 
far  at  present. 

With  the  introduction  of  an  additional  hypothesis  governing  crack  growth  the  present  model 
can  be  used  as  the  basis  for  a  theory  of  tearing.  That  something  about  the  shape  of  the  crack  tip 
remains  invariant  as  the  crack  grows  has  been  suggested  as  such  an  hypothesis  [9,  10].  An 
adaptation  of  that  hypothesis  to  the  circumferential  crack  problem  will  be  the  subject  of  a 
subsequent  paper. 
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